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SUMMARY 



In order to determine the effect of a non-linear 

spring characteristic on the vibrations of a torsional 

<• 

system, investigations were made to determine the response 
carves of a given system when the drive was made through 
a non-linear coupling and when made through a comparable 
shaft* The results showed the maximum amplitude for the 
former case to be about one-fourth of that for the latter. 
However* when the difference In damping was acounted for, 
the non-linear coupling showed a reduction In amplitude 
of vibration to about one-half that of the comparable 
linear system. An analytical solution including damping 
is given, and suggestions are made for further study. 
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INTRODUCTION 



While couplings having non-linear torque-def lec- 
tion characteristics have boon used for some time In In- 
stallations where It is desired to reduce torsional vi- 
brations* and while the accomplishment of this mission 
Is generally conceded, there la no data available for 
evaluating their effectiveness In this regard. It la de- 
sirable, when considering an Installation of this type, 
to know not only if an improvement will be realized, but 
also how much improvement. It was, therefore, with this 
point in mind that the subject investigation was under- 
taken. 
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Pig. 1 




Pig. 2 



3 



A — ammeter (load amps.) 

B — pedestal bearing 
a - b — hitworth drive 

b — driving disc (referred to coupling) , graduated 
0-300 degrees 

c - c — non-linear coupling 

<1 — driven disc (referred to coupling), graduated 
0-360 degrees 

G — - generator (d.c., 115 v,, l/20 H,P. at 1725 R.P.Sf.) 
L — lampbank load 

M — motor (d.c., 115 v., 2,43 amps,, l/4 H.P. at 
1725 R.P.M.) 



P — - breaker points, attached to graduated dial 
R — rheostats for speed control 
r — rheostat for control of generator field 
Ra — mounting rail 
S — stroboscope 
T — torsiograph 

Y|.— voltmeter (generator field voltage) 

V t — voltmeter (generator terminal voltage) 

X — ?, strobotac” (not shown) 



Powori 6 v. d.c,— for atroboscopo and torsiograph 

115 v. d.c. — for motor drive and generator field 
115 v. 60 cycle a.c.— for strobotac 
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Fig. 3 




Fig. 4 
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The Whitworth drive is an inversion of the 3llder 
crank mechanism designed to impress an oscillatory torque 
on disc ”b". The offset between the shaft to "a" and the 
shaft from "b" was kept constant at l/8" throughout. 

The coupling used is shown in detail in Figures 3 
and 4. It may be seen that the segments holding the spring 
strips are cut with first a plane portion and then a cylin- 
drical portion. The radius of curvature of this latter part 
is 6 inches. The non-linearity is produced by the wrapping 
of the spring strips upon the cylindrical segments thus: 



This shortens the unsupported length of the strip and so 
stiffens the coupling. The spring strips have a free slid- 
ing fit In the straight portion of the segments in order to 
allow for the axial motion necessary when the coupling is 
deflected. Clock-spring stool .025 inches thick and 3/8 
inches wide was used for the spring strips. The "coupling 
lengths" referred to In the discussion to follow are the 
unsupported lengths of the spring strips when not deflect- 
ed. Thi3 length wa3 varied for the different runs by slid- 
ing the assembly "G", "d", and one coupling flange ”c" to- 
ward the motor end of the set-up. The ends of the strips 
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wore then cut off to allow only the alight projection beyond 
the coupling flanges shown in Figure 2. Parts "M* , "B" , end 
"G* were bolted rigidly to the mounting rail M Ra M . For the 
runs made with ordinary linear shafts only the spring 
strips wore removed from the set-up as shown In Figure 1. 

The coupling flanges M c" were allowed to remain In place 
and the shaft connected between them by means of chucks. 

Two Intermediate bearings were used to support the shafts 
tested. 

The shaft of motor ”W M projects through the out- 
board bearing to within a short distance of disc "P”. A cam 
on this shaft operates a set of breaker points mounted on 
M P”, breaking the circuit through these points once each 
revolution. The primary of the stroboscope is in series 
with this circuit, causing the stroboscope to flash each 
time the circuit is broken. By rotation of disc "P M to the 
desired position this flash may be made to occur at any 
position of the shaft during a revolution* A knurled lock- 
ing nut is used to lock disc n P* in the position desired. 

The effect of this mechanism Is to allow simultaneous read- 
ings to be taken of the position of discs H b” and ”d w at 
various points throughout a revolution of fl b” , these points 
being determined by the setting of dial M P M which is gradu- 
ated. Thus a scan may be made of these instantaneous posi- 
tions throughout a revolution. A circuit from the breaker 



points also operates the dead center timing pen of the 
toraiograph. 
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The toraiograph is driven by a belt from the shaft- 
ing between "b" and "B”. In order to reduce slippage this 
section of shafting vms wound with friction tape. In this 
condition the shaft diameter was ll/l6 Inches. ¥<hen run- 
ning the linear shafts it was found that the amplitudes of 
vibration were too large to be rocorded by the toraiograph 
as installed, even at its lowest amplification. This was 
remedied by driving the toraiograph from the test shaft It- 
self, which provided the necessary smaller driving diameter. 
This shaft, taped, measured 7/32 inches. 

The generator field was separately excited, the ex- 
citation being maintained at 100 volts. The lampbank load 
was not changed during the investigation; hence the load In- 
creased with the speed of operation. The lamps of the lamp- 
bank were painted black in order that they would omit no 
light . 

The strobotac was mounted in the position Indicated 
by M X" In Figure 1, so as to throw Its light on the spring 
strips. hen running the shafts, the strobotac was brought 
to the near side (In Figure 1) of the mounting rail and Its 
light directed on that face of disc w d” which faces the 
coupling flanges. The strobotac was used for speed measure- 
ment only , and should not be confused with the stroboscope 
which was used In obtaining dial readings as described above. 
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For all runs Involving the use of the stroboscope 
or the strobotac, the laboratory was made as dark as pos- 
sible In order to facilitate the taking of readings. 



PROCEDURE 



With the apparatus sot up as shown in Figure 1, a 
preliminary run was made to determine the speeds at which 
it would be desirable to take readings. This was done by 
increasing the motor speed by small increments, each time 
"stopping" the rotational motion of the coupling with the 
strobotac. In this state the amplitude of the vibratory 
motion in the coupling could be observed by eye. Speeds 
were thus chosen which gave desired increments of ampli- 
tude change. The critical speeds were also determined 
during this run. 

The strobotac was then sot at the lowest speed at 
which an amplitude determination was desired. The motor 
speed was adjusted, by means of the motor rheostats, until 
rotational motion of the coupling was "stopped" by the stro- 
botac light. Since the speed of the motor had a tendency to 
vary, it was watched constantly and readings were taken only 
when the speed was steady and at the desired value. The ro- 
tation of the motor shaft caused the stroboscope to flash 
once each revolution, as previously explained. By the 
light of those flashes the dials on the driving and driven 
discs (designated as "b" and "d" in Figure 1) were read 
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and recorded. The corresponding setting of dial n P” was 
also recorded. Dial "P n was then rotated in increments of 
15 degrees, thu3 causing the stroboscope to flash at suc- 
cessive positions during the cycle. Readings of dials "b” 
and "d" were taken for each position. When a complete cycle 
of readings had been obtained the voltmeter and ammeter read 
ings were recorded and a record of the motion of disc "b" 
taken on the torslograph. The speed was then increased to 
the next desired value and the procedure repeated. Readings 
were thus taken at 3peeds ranging as close to the critical 
speed as possible. Readings at the critical speed were not 
possible due to the difficulty of maintaining constant speed 
After passing the critical, the speed was Increased to the 
maximum value at which readings were desired, and then re- 
duced in steps, taking readings at each step as before. 

This was continued until the point of Instability at de- 
creasing speed was reached. 

The motor was then stopped and a reading taken of 
the relative position of dials n b” and ”d lt with no twist 
in the coupling. This served as a reference reading in 
determining the actual position of one disc, relative to 
the other, from the readings previously taken. 

Following this a bar was attached to the shaft be- 
tween w B n and "c" and the end of the generator shaft was 
rigidly clamped to prevent turning, by means of a pedestal 
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secured to the mounting rail. Using a spring scalo known 
increments of torque wore applied to the shaft and the re- 
sulting angular deflection read from dial "b". This sas 
done for increasing and decreasing values of torque to ob- 
tain data for the hysteresis loop of the unit. This com- 
pleted the readings taken for a given coupling length. 

The entire procedure was repeated for each of the 
coupling lengths. 

The procedure for each of the shafts (i.e., "linoar 
couplings”) tested was identical, with the exception that 
the hysteresis loop was taken before the run in order to 
determine that the shaft stiffness was of the desired mag- 
nitude. This stiffness was also checked after the run to 
assure that no change had occurred due to the stress of the 
largo deflections encountered. 

Due to the small output of the generator and the 
difficulty of determining its efficiency at this light load, 
a series of runs was made with the motor and generator 
coupled through a short, stiff shaft. The input to the 
motor was determined with no load and no field on the 
generator, and again determined with the load and field 
adjusted to that used in making the coupling runs. The 
difference w&3 assumed to be the load input to the gener- 
ator. This difference was determined for various speeds 
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and the results plotted (see Figure 18). **hll© the accuracy 
of this procedure is not good. It is considered satisfactory 
inasmuch as it contributes only a small fraction of exciting 
torque applied to the driving disc. 
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THEORY 



1. Notation 



a 



e 

I 



T 58 

i E 



component of steady state amplitude, see defini- 
tion of a and y. 

viscous damping factor, pound Inch seconds per 

radian. 

component of steady state amplitude, see defini- 
tion of s and y. 

the eccentricity between discs *'a tf and M b” . 
moment of inertia of mass or combination of 
masses Indicated by subscript, see Figured, 
mass moment of inertia of equivalent system 



k - 



ki a 

k 2 = 



torsional atiffnoss of shaft, pound inches per 
radian. 

initial slope of non-linear torque-deflection curve, 
coefficient of the cubed term in the expression 
for the non-linear torque-deflection curve. 
k 0 



3/4 



n 



N 



Jb 

21 



E 



30 



rotational speed in revolutions per minute * ~ *f 
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N 



n 



p m 



natural frequency expressed In revolutions per 

minute * ~ 

IT 

r 

W 



r * the crank radius of the Fhitworth linkage. 

R * frequency ratio “ -i * 2— 

*n N n 



S * y^ * 

T ” f T ♦ q» 

*a L 

Tb* the coefficient of the sinusoidal component of 
driving torque applied to disc "b". 

Tb'“ total driving torque applied to disc "b”. 

Tjv 3 * damping torque. 

Ii Ip 

Tg* equivalent torque “ — = 

X 1 * *2 

Tj“ load torque. 



Tj® Inertia torque. 

*r T i ♦ T b 



y * angle of twist in radians. 



l E 



st kf 
ft m “amplification factor” 






dlmensionsless damping coefficient 



/(1-R^) S ♦ 'k s 

b Zn 



I U) w 

n n 
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** angle of twist in degrees. 

*f" circular frequency of forced vibration. 
w n * natural frequency *7/ ^1 
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2. Derivation of amplitude equation 



In order to reduce the complexity of the analysis 
the system must first be simplified. The steps of this 
simplification may be indicated as follows: 








E 







Figure 5 



The letters indicating the original masses of the system 
above are the same as those employed to designate these 
parts In Figures 1 and 2. 



Since the masses w M n and ”a" aro rigidly connected. 



I 



A 






* i. 



likewise. 



h - h - *c 
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and I 2 - I c ♦ I d * I G 



Equating the kinetic energies of system 1 and Its component 
parto: 

1/2 I A ♦ 1/2 I B «§ * 1/2 1 1 «av. 

av* av* 



But a) 



av< 



% 



<*> 



av, 



av< 



. . - I A ♦ ! b 



Y<e now have a two mass system with moments of inertia 1^ 
and I g , connected by a spring of stiffness k. 

Let a periodic torque of magnitude be applied 
to mass 1. Then It may be shown (Reference 1, page 33) that 
the equation of motion of the two mass system Is: 



I 1 I 2 y 
L 1 ' A 2 



+ by * ky - I 2 T 1 f (»,t) 

+ T 

1 x 2 



It + Io I_ + I 



Where y * the relative motion between discs M b M and "d M 
(viz., the angle of twist). 

This equation Is Identical with the equation of absolute 
motion of the system E If 



x l x 2 



X 1 * X 2 



and Ttr* X 2~l 

'• v — y f" ~ 

x i H 
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Therefore we may write as the equation of motion of the 
equivalent system: 

I E y + by ♦ ky * T £ f(«> f t) (1) 



T_ f(w t) may be written as a Fourier series, the motion 
due to each component evaluated, and the net resultant mo- 
tion obtained by superposition. From this point on we there 
fore will consider only the fundamental component T„ sin <*>t 
which satisfactorily describes this torque for the case 
under consideration. 

In the case of a linear system (k * constant), wo 
have a linear differential equation, which may be solved 
for y. giving the well known expression for maximum ampli- 
tude: 



TvA 






max 



/(l-R 2 ) 2 ♦ sV 



- (3 (t e A) 



However, for the non-linear syaten [k * f(y)J no 
such simple general solution exists. The following deriva- 
tion Is based on an article by Appleton (Reference 2). 

First, it la assumed that the torque deflection 
curve of the torsional spring can be expressed by an equa- 
tion of the form: 

T “ k x y * kgy 3 

The choice of this form of expression for the spring charac- 
teristic is arbitrary, and Is chosen primarily for conveni- 
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lence. However, results of this Investigation show It to 
give satisfactory agreement with the actual curves obtained. 

Substituting this expression for our spring torque 
Into equation (1) we have: 

I K y ♦ by ♦ ^y ♦ k 2 y° - T E sin # f t, (2) 



or 



2 . 2 , 



I g y + by ♦ k^y (1 ♦ jr- y ) * T g sin . 

1 

Dividing by I £ , we obtain: 

b k l kg 

y + — y + ij“y(i ♦ ej y ) * sin « f t 



L E 



or j 



Y + 



ko <> 



£ny ♦ w n y (1 ♦ y") - w n y st sln li f t » 



where 2n = 



n 



h 



and y 



st 



T i 

*1 



Neglecting the transient term, let the trial solution be: 
y * a sin « ,t + c cos u>,t (4) 

then y = a u> .cos - c m £ sin w t (5) 
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and j m - a sin W^t - c cos u»j,t (6) 

3 

In considering the non-linear term involving y , 
we shall, as a first approximation, neglect the influence 
of the harmonics and retain only those terms involving 

Thus .cubing equation ( 4 ), reducing terms to first 
power functions of w t. and discarding all functions of 
t, 3 u>.»t, we obtain; 

y 3 * S/4 a 3 sin ♦ 3/4 a‘c cos «tfft + 3/4 ac sin fc^t 
♦ 3/4 C COS 

Substituting f t y and y 3 from equations (6), (0), and (7) 
into equation (3) and equating the coefficients of the 
sine and cosine terras respectively, we have; 



-aw - 2ncct>. + w'a + 3/4 a ^2 (a^ ♦ c ) « y g< . 

k l 

-cwf + 2naw f ♦ #n c + 3/4 c k 2 (a' ♦ c ) • 0 



( 8 ) 



Rewriting: 

a - 2nm c + 3/4 

c (w^ - tt^) ♦ 2n»|.a + 3/4 



k 2 a (a' ♦ c 2 ) - wf y at 
k l 

^ c (a 2 ♦ c S ) « 0 



( 9 ) 
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U> f 

Let R ~ , 

w n 



X - , 

k 2 

K - 3/4 ^ . 

Substituting these relationship® Into (9), wo obtain: 

a(l-R 2 ) - J'Rc ♦ aX(a* + c 2 ) * y Qt 

, 2, v .2 2 

c ( 1-R ) + Vila * cK( a ♦ c ) - 0 

Rearranging: 

[ 1 - K 2 ♦ K (a 2 ♦ c 2 ) ] a - (fRc - f at 

{ 10 ) 

f 1 - R 2 + K (a + c c ) ] c ♦ ^Ra * 0 

2 2 

Let S «* a*' + c 3 the square of the amplitude for steady 
state condition*. 

Substituting and operating on (10) : 



(1 - R 2 ♦ KS) 2 a 2 - 2(1 - R 2 ♦ KS)a ^Rc ♦ tf 2 R 2 c 2 - y 2 t 
(1 - R 2 ♦ KS) 2 C 2 ♦ 2(1 - R 2 ♦ KS)c ^fRa ♦ fr 2 R 2 a 2 - 0 



Adding: P „ 

(1 - R ♦ KS; 5 ♦ frTT* - j% t 



( 11 ) 



This is a cubic equation in S, the real roots of which will 
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give values of the possible steady state amplitudes. The 

criteria for stability of these amplitudes will not bo dis 

cussed here since the regions of stability and instability 

of the non-linear response curve are well known and are 

covered in most vibration texts. 

Since solution of a cubic equation is unhandy, S 

nay be made the independent variable and the equation ox- 

2 

pressed as a quadratic in R . 

Rewriting equation (11): 



[ (1 + KS) - R 2 ] 2 + - y 2 

S V; 



S 



(1 ♦ KS) 2 - 2(1 • KS)R 2 + R 4 ♦ <f 2 R 2 - y 2 t/ , - 0 



R 4 - 2(1 ♦ KS - * 2 / 



% >k 2 - y»t/ g - (i + 



KS) 2 - 0 



whence 



R 2 - KS+d-y 2 ^) 1 / K 2 3 2 + 2(1-^) KS+(l-^4 f +y 2 t - 1- 2KS ~ K^S 2 

'S 




Equation (12) is the expression used in the compu- 



tation of the non-linear coupling response curves. It lends 
Itself readily to a computing form. 
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3. Other theoretical considerations 



(a) Evaluation of 

Prom the preceding derivation; 

„ T l I g 

E *1 * I 2 

The net torque exerted by the motor must be equal 

to the resisting torque of the generator plus the damping 

torque, T and T respectively. This net torque is a con- 
L D 

stant at disc "a”, but In being transmitted through the 
linkage of the Whitworth drive acquires an oscillatory com- 
ponent. Prom the mechanics of the system it may be deter- 
mined that 

T ’ - T a P 2 jS 2P ' C ° 8 ‘fV - * 
h a p - p cos (wt) 

where p - §• 



r *» the distance of the driving pin of the linkage 
from the center of disc "a", equal to 1.5 inches 
for this installation. 

e » the eccentricity between discs n a" and ”b n , equal 
to 0.125 inches for all runs made. 

V* t l + T D 

“ total driving torque applied to disc ,! b”. 
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Evaluating the constants and expressing as a series 

wo have: 

Tjj’ * T & ( 1*00- .083 cos tat - ,003 C 03 2a>t - .0001 cos 3*>t 

♦ ..) 

Discarding the harmonics above the first as negli- 
gible, there remains, superimposed on the constant torque, 
a simple sinusoidal torque of maximum amplitude: 

T b - .083 {T l ♦ T d ) 

Since vibrations of the masses ”B n may be trans- 
mitted through the Whitworth drive to masses ”A W , the latter 
will be caused to oscillate also. These oscillations result 
in an inertia torque being Impressed on disc "b". Desig- 
nating this torque as T * we may write: 

T t * T o( 

l I a A^A 

Again by the mechanics of the Ahitworth drive, sub- 
stituting the constants used above, we have: 

A 13 (.083 sin *t • .007 sin 2#t • .0003 sin 3 fit 

♦ ) 

Neglecting the harmonics higher than the first and 
substituting above we obtain: 

Tj* » .083 I A sin Wt 



. 














. . . 
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Or, designating the maximum amplitude of the sinus- 
oidal torque as Tj, and noting that v» the same as the 
speed of rotation of the coupling (via., ) , wc have: 

T 2 * .083 W f I A 

The oscillatory component of the torque applied 
to system ”l w will be the sum of the oscillatory component 
of the driving torque plus the inertia torque. This will 
be a sinusoidal torque of maximum amplitude: 

T x - Tj ♦ T b 

Substituting in the expression for Tj , we have: 

s - [- 085 •»* h • -° 83 < t l ♦ V] r-?V 

X 2 

- [* 083 **■ * V * - 083 (T l * T d )] 

Tj Is evaluated from a speed vs. load-torque curve. 

T d is evaluated from the hysteresis loop for the par- 
ticular set-up being considered. 

(b) Evaluation of the damping factor, ^ . 

The damping factors used in the computations for 
both the shafts and the non-linear couplings were so deter- 



. 



















. 










. 

. 

■ 
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mined as to make the maximum amplitude of the theoretical 
curves at resonance the same as that observed. Although 
this is apparently a subterfuge, its use will be defended 
in the discussion to follow. 

For the linear system at resonance, i and 
may be determined by computing backwards from a given maxi- 
mum amplitude. 

For the non-linear system 0 may be determined as 

follows: 

Differentiating equation (11) with respect to R, we obtain: 

[(1 - R 2 ) ♦ XS] 2 H ♦ 2S [(1- R 2 ) ♦ KS] (-2R ♦ Kg|) 

+ K 2 R 2 -4^ + 2 tf 2 R3 * 0 

dR 

[(1 - R 2 ) 2 ♦ 2 KS (1 - R 2 ) ♦ K 2 S 2 * 2 KS ( 1 - R 2 ) 

+ 2 tV + *VJ - 4 RS [Cl - R 2 ) ♦ K* • If ] 

or: 

2 KS (2-2 R 2 + 2 XS - ’) 

* 5—5 - « 

dR (1 - R ) + fR ♦ 4 FS (1 - R^) + 3 K S 

The maximum point on the resonance curve occurs when the 
sldpe la zero. Since R and S are not zero, it is neces- 
sary that: 



0 
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2 - 2 R 2 ♦ 2 KS - ft - 0 

or 

1 - R 2 ♦ KS - $ 2 / g - 0 

Substituting for R 2 from equation (12): 

i - m - i • ^/ 2 ; ^ KS+ < i- ) 5 -i ♦Ks-jP/g • 0 

/ s 

whence 

y?t/ - * 2 KS ♦ (1 - 5%) 2 - 1 - 0 

y 2 t/ - fr S KS ♦ 1 - & ♦ - 1 “ 0 

'S 

S 4 /* - (KS ♦ 1) ft" ♦ y 2 y - 0 

2 ( KS ♦ 1 t XkS +1) 2 - y 2 . 

’s 



> 



( 13 ) 



DISCUSSION 



Since the comparison of the non-linear with the 
linear system wa3 made on the basis of an ‘'equivalent shaft” 
it is necessary first to define this term. This definition 
follows: "A shaft is considered equivalent to a given 

coupling when the slope of its torque-deflection curve is 
the same as the initial slope of the torque-deflection curve 
of the given coupling.” On this basis the Initial flexi- 
bility of a coupling and its equivalent shaft is the same 
(l.e., k of the shaft equals of the coupling). Thus the 
two systems are brought together at an initial point. De- 
parture of the behavior of one system from that of the 
other as operation occurs removed from this point (viz., 
amplitude of vibrations become greater than zero) must be 
due to whatever difference exists in the spring character- 
istics and the damping of the two systems, since the moments 
of inertia remain essentially the same. 

As the purpose of this investigation was to deter- 
mine the effect of the non-linear spring characteristic, 
it would have been desirable to make the damping of the two 
systems identical* However, no attempt was made to do this, 
and it will therefore be necessary to estimate the effect 
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of the difference in damping In order to draw a conclusion 
as to the effect due solely to the non-linearity of the 
spring. It is doubtful that a coupling could have been de- 
signed such that the damping in the coupling system and in 
the equivalent shaft system would have been identical. In 
neither case was the damping purely viscous; rather It was 
of some complex type, evaluation of which would have been 
quite difficult. Following accepted practice, the damping 
was expressed as an equivalent viscous damping, in amount 
equal to that necessary to give the observed maximum ampli- 
tude at resonance. The method of making this evaluation 
for each system has bean explained under "theoretical con- 
siderations. w It is believed that the major cause of dis- 
crepancy between theoretical and observed values of ampli- 
tude Is due to this approximate means of evaluating the damp- 
ing, and it is therefore suggested that if future Investiga- 
tions are undertaken, they be directed primarily toward a 
more accurate solution of this problem. 

Accepting the damping determined os described above, 
it was still necessary to arrive at some basis of comparison 
which would present the true effect of the non-linear spring 
characteristic alone. Since the damping in the coupling was 
much greater than that In the shaft. It Is obvious that a 
comparison neglecting the difference in damping would present 
the former in too favorable a light. Therefore a theoretical 
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maximum amplitude was computed using the valueo of damping 
and torquo at resonance of the equivalent shaft. The maxi- 
mum amplitude thus obtained was almost double that observed, 
but still about half the maximum amplitude of the equivalent 
shaft, as may be seen In Figures G to 9 inclusive. 

Damping also complicated the evaluation of the 
torque-deflection curve of the non-linear system. Since 
the deflection was changed very slowly In obtaining the 
curve, the only damping present was that Independent of the 
speed, viz.. Coulomb damping. The result of its presence, 
however, was to give a hysteresis loop rather than a single 
curve. Inspection of the curve indicates that the friction 
was greater with decreasing torque than with increasing 
torque. Thus we have tho spring characteristic confused 
with part of the damping. The true spring characteristic 
must lie somewhere within the loop. A characteristic evalu- 
ated to pass through the maximum point of the hysteresis 
loop was found to be unsatisfactory since, at the larger 
angles of twist, conditions are such as to make the damping 
effdet overshadow the spring effect on the curve. Similar- 
ly, the mean line between tho two boundaries of the loop was 
found to give poor results because it was unduly lowered by 
the increased damping effect of the descending branch of the 
curve. By trial It was found that a curve approximately 
parallel to the upper branch of the hysteresis loop, passing 
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through the origin, gave the most uniformly satisfactory 
results. Since the initial slope of the approximate curve 
represented by the cubic equation is fixed (it is equal to 
the initial slope of the upper branch of the hysteresis 
loop), it can be forced to coincide with the actual curve 
at only one point other than the origin. This coincidence 
was made approximately at the mid-point of the actual torque- 
deflection curve. Since the actual curve itself is at best 
only an estimate as shown above, a small departure from it 
by the approximate curve is not considered serious. 

Referring to equation (12), it should be noted that 
T k Is a function of (or of B), which Is In turn a func- 
tion of R. The equation is thu3 not In proper fom for 
evaluation as a quadratic. However, a value of H may be 
estimated, from which may be computed. If the value of 
R thus determined gives an H which differs by more than 10 
per cent from the M estimated, the computed N may be used 
for the next estimate and the value of R recomputed. If 
there Is no other means of estimating R for the first ap- 
proximation, the value of N n (which is known for any par- 
ticular system) may be used for all but the lowest ampli- 
tudes. If the first estimate of H is considerably lower 
than the actual value, the quantity under the radical be- 
comes negative, giving a complex number for R 2 . In general, 
knowing tbe value of N n , and having a knowledge of the gen- 
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eral shape of non-linear response curves, a sufficiently 
satisfactory estimate can bo made so that not more than two 
approximations are necessary. A 10 per cent variation in 
the estimated value of H will give approximately a I per 
cent variation in th© computed value of N. Prom this it is 
apparent that the result is not being forced by the esti- 
mate. If equation (// ) had been solved as a cubic in S It 
would not have been necessary to make the above approxima- 
tions since in that case R would be the independent variable. 
It is considered, however, that the above method lends it- 
self more readily to mass calculation (i.o., a computing 
form) • 

The observed amplitudes indicated on the response 
curves. Figures <3 to 9 inclusive, were obtained directly 
from the readings of dials n b !> and "d" and represent th© 
angle of twist from the mean position, viz., half the sum 
of the maximum angles of twist in each direction. 

It is interesting to note that in the cases whore 
the non-linearity is most pronounced, the overlap of the 
upper and lower branches of the response curves is suffi- 
cient to allow advantage? to be taken of the effect, as fol- 
lows: If it is desired to run at a speed approaching the 
critical, increase the speed first until the critical Is 
passed, then reduce the speed to that desired. Thus the 
machine may be made to operate on the lower branch of the 
curve where the amplitude is much smaller. 








. 






CONCLUSIONS 



The results of this investigation indicate that a 
non-linear spring characteristic in a torsional system is 
effective in reducing the maximum amplitude of vibrations. 

An exact determination of the magnitude of this reduction 
was complicated by the existence of complex damping in the 
systems compared. However, a reasonable approximation seems 
to be that for systems of tho type investigated the non- 
linearity will cause a reduction of about 50 per cent in 
tho maximum amplitude from that of the comparable linear 
system. 

It was also demonstrated that in a non-linear sys- 
tem advantage may be taken of the overlap of tho upper and 
lower branches of the response curves to further reduce the 
maximum amplitude which must be accepted when operating In 
the critical zone. For the cases considered this additional 
reduction amounted to about 80 per cent of the maximum ampli- 
tude Indicated by the non-linear response curve. Referring 
to Figures 6 to 9 Inclusive, it Is seen that the response 
curve of the non-linear system shows appreciable amplitude 
over a wider range of frequencies than does that of the 
linear system. However, tho great reduction effected In 
peak amplitude is considered to more than compensate for 
this disadvantage in most installations. 
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APPENDIX I 



Sample Calculations 

In the following we shall consider first, some gen- 
eral calculations used for all computations. Following 
this, sample calculations for the 4-l/o inch coupling are 
presented, together with those for Its equivalent shaft. 
Finally, computations are made for the theoretical maxi- 
mum amplitude considering the shaft torque and damping to 
be applied to the coupling system. 



General Calculations 



I - I ♦ 3 

1 A S3 



I ♦ T ♦ T ■* T 
"a x b A c 



- .01709 ♦ .02650 + .02354 «■ .01085 



** .078 slug In.' 



' x c + *d + *0 



“ .01085 ♦ .06594 ♦ .00125 



* .078 slug in. 
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. ^Vu-rJ I 

































.078 x .078 _ _ ? 

*» . 03 9 slug ln . 

*2 T 1 . *078 T x T i 

I x ♦ Ig * .078 ♦ .078 " 2 

Calculations for the 4-1/3 Inch coupling 

Evaluation of the approximate torque-deflection curve. 

From Figure 11, the initial slope of the torque- 
deflection curve is: 

k x - .53# In./o x 57.3 « 30.4# in. /rad. 

and T * 17. 6£ in. when 0 * 20°. 

20 

y “ ■ -- * .349 radians 

Of* 4 s 

kp 3 

? * k x (y + y ) 
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Determination of the natural frequency. 



*h - ni 

X K 



ySTTi T"~ =* 28 rad. /see. 



H 



n 



30 

7T 



c V 



n 



268 r.p.m. 



Evaluation of the damping* 

The maximum amplitude observed for this coupling 
was 29*5° at H * 387. 

(Vx, * JL u » 13 40.5 rad*/sec. 

i 30 30 



T* 3 *083 [ t* f (I K ♦ I a ) ♦ (T t . ♦ T n ) ] 



D' 



T + T 
1 A 2 



- .083 (TO .5 (.01709 + .02650) ♦ (2.50 ♦ 1.20) ]x 



“ 3. IE# in. 



y 8 t 

2 

y.t 

ymax 



, T K * 5.15 .. 

*s — - ms — , , - , 

k x 30.4 
* .0107 



29.5 

57.3 



.514 



.1036 



S 

max 



2 

y 3 t 



s max 



^ O 

y " = 70X4 * .264 

J max 

* .0107 » .0404 

TOOT* 



CO 



■ 

. . . 









. 












. 



■ 









- 






■ 













: 
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ft 2 - 2 (KS ♦ 1 t V (KS + 1)2 _ y| t ^ ) (13) 

Substituting: 

^ 2 =2 (4.1 x .264 + 1 + y(4.1 x .264 + l) 2 - .0404 ) 
ft 2 * .0196 
)f = .1400 

Evaluating the response curve: 

r 2 = ks + (i - *%) i V ylt/ s - + (l - &%)* " 1 U2> 

Substituting: 

R 2 - 4. IS + .990 ♦ Vy - .08S - .020 

923S 

o 

VJe may now substitute values for S and solve for R . 

For convenience, take the amplitude y = .1 radian (ft - 5°.73). 
Then S = y 2 = .01 

To evaluate T^, it is necessary to estimate a value of N. 

A substantiation of this procedure is given in the discussion. 

For y * .1 , for the upper branch of the response curve 

N - 220 RPM 
est 

For the lower branch of the response curve 



370 RPM 



40 



Proceeding for the upper branch of the response curve: 

o Io 



T e - .083 <I„ + !„)♦ (T t ♦ T d ) ] 



T ♦ T 
X 1 X 2 



» .083 [ (£|~) (.0436) + 2.76 ] xi 



30 



- 1.09 



Tj and were obtained from Pigs. 18 and 11 respectively. 



o 

Substituting in the equation for R and using the negative 
sign before the radical since the lower of the two values 
for R 2 is desired: 

R 2 » 4.1 x .01 + .990 — - .08 x .01 - .020 

v 923 x .01 

- 1.031 - .3289 

» .7021 

R - /.'TOS'I - .837 

K - R}1 = .837 X 268 » 225 RPM 

n 

Proceeding similarly for the lower branch of the response > 

curve, and using the positive sign before the radical, we 
obtain: 

R - 374 

Thus ttoo points on the response curve have been obtained. 




■ 







- 



* 



. 
























. 

. 




























41 



viz . : 



« 5?75, N * 225 



5?73, M - 374 

Following the same procedure sufficient points may he ob- 
tained to determine the theoretical response curve for the 
4-l/3 inch coupling given in Fig. 7. 

Computations for the equivalent shaft. 

Using calculations given above for the coupling: 



k « k-^ * 30.4 ’* n / rad 



U> n * 27.8 rad/sec 



N 3 208 
n 



At resonance 



T j, ** .083 



[<^M .0436) + 3.08 ] X 1/: 



2 * 1.54 



0 obs. * 100° 



y 31 * 1.747 rad 

v> Y • u 



R - 1.0 








. 

. 



































. 
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1.747 



1.54 
30.4 K 



X * .0292 



Values for the response curve may now be computed: 

Determine the amplitude at R * 250: 

T« “ 1.35 , R » 552 - .933 , R 2 - .870 

* 268 

T Vk 

^■‘ a X V" " """ff '- g r £>" '"& — 

Mi - ) 2 + r 2 ir 

1,,55 /30.4 

g 2 

(1 - .87) ♦ .0202 x .870 



“ .334 rads. 

« .334 x 57.3 * 19?1 

This gives one point on the linear response curve of Fig. 7. 
Other points may be determined similarly. 

Computations for tho theoretical maximum amplitude 
of the non-linear system when the torque and damping at 
resonance of its equivalent linear system is applied to it. 
From previous computations - 
T e - 1.54 

« .0292 



Under ” theoretical considerations" it was shown that 




. 



at maximum amplitude of the non-linear response curve 



S 4 / 4 - (KS +1) tf 2 * yg t / s - 0 

Substituting and solving for S - 

- . 4 2 

.0252 - (4. IS ♦ 1) x .0292 

— r~ 

S - .728 

y *= VS - * .853 radians 



+ 1.54 - 0 

dS'Ss 






^ • •• f k -fc 
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Data 



Momenta 


of 


inertia: 






** 


2S 


.01709 


slug 


In 


*0 


a 


.00125 


If 


it 


X a 


as 


.02650 


If 


w 


h 


* 


.02354 


ft 


ft 


*c 


as 


.01085 


If 


n 


J d 


<25 


.06594 


n 


if 


"Equivalent 


" shafts 


• 

# 





All "equivalent” shafts were of mild 3teel, 3/l6 
Inches diameter* The lengths wore as follows: 
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Inch 


shaft 


equivalent 


to 


4 


Inch coupling 


45-1/2 


It 


n 


H 


t* 


4-1/3 


ft « 


50-5/8 


M 


n 


U 


« 


4-2/3 


n it 


60-7/8 


n 


n 


ft 


II 


5 


« ft 
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T I.11AT IO* OF CtrfEOJSKD B'TA AID CWS^iBIXC Cfj&rpUTFB VLUES 



4 inch Coupling 


4-1/3 inch Coupling 




ob*#rrod 


Computed 




Observed 


Computed 


Frecueney 


Amplitude 


Amplitude 


Frequency 


Asplitude 


Amplitude 




0-dogree* 




S-tUTf 


’0 -decrees 


i 

£ 


250ft 


4.5 


5.8 


200ft 


3.5 


3.5 


300 f 


11*5 


11.5 


250ft 


0.0 


9.5 


340ft 


17.5 


16.3 


230ft 


14.0 


14.5 


350ft 


13.5 


17.5 


300ft 


19.5 


17.7 


4094 


22.0 


23.5 


320*. 


23.5 


20.7 


425 ft 


25.0 


26.0 


340 ft 


25.5 


23.7 


430’ 


28.0 


2G*3 


3J0ft 


2C.5 


28.5 


330D 


10.0 


9.0 


370ft 


28.0 


28.0 


400D 


7.0 


7.0 


330" 


29.5 


23.1 


4S0D 


4.6 


5.0 


333ft 


29.5 


29.4 








S40D 


12.5 


9.5 








35013 


8.5 


0.0 








400D 


6.0 


5.0 


| 


39 inch Shsj 


“t 


4. 


r*l/2 Inch S 






Observed 


Computed 




Observed 


Computed 


Frequency 


Amplitude 


Acaplitudo 


frequency 


Amplitude 


Amplitude 


""MS ! 


0-degr*®* _ 


^-degrees 


iM»Si 


^'-degrees 




250 


8.0 


7.7 


200 


3.5 


3.5 


235 


13.0 


12.5 


240 


12.0 


12.0 


275 


27.5 


18.5 


250 


20.5 


19.7 


283 


5S.5 


37.5 


200 


47.5 


55.0 


305 


16.5 


32.5 


2S5 


18*0 


13.3 


325 


11.0 


15.5 


300 


12.6 


13.3 


370 1 


7.5 ' 


7.5 


400 1 


5.0 


5.0 


400 


5.0 


5.7 








450 


4.0 


4.0 









lota: "A* indicates operation on uppor branch of non— linear response curve. 



*B* indicates operation on loser branch of norr-l irtenr response curve 
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TA3UL*TZ < M QT COHFMSW D'fA S3 vO**' - . ODIIU vC’^lfl® ?*U«< 



4 - 2/3 inch Coupling 


5 inch templing 


frequency 


Obeervod 


computed 


ne-uenoy * 


| tbeervod 


Computed 


>«plltud» 


"jsplitbde 


rplitucle 


t-plltude 






0-d®cwN»* 


h 1 


| P'Klil 




200' 


3.5 


4.6 


200 


3.5 


C.4 


230' 


7.0 


0.2 


225' 


8.0 


1 .5 


2S0‘ 


14.5 


14.5 


2501 

275* 


17.0 


14.7 


260 


24.0 


13.8 


19.S 


10.0 


31C6 


26.0 


23.2 


3CD» 


22.5 


23.3 


525‘ 


2C.5 


25. C> 


510* 


25.5 


24.5 


336“ 


26*0 


27.0 


520* 


23.5 


25.5 


31Sfi 


».e 


0.0 


90PD 


0.0 


8.2 


ssor, 


*.5 


5.5 


340£ 


0.0 


5.8 




4*0 


4.0 


3050 

1 - ! ! t | . . -1 


8tC 


4.8 



fO*/k inch 



:h*rt 



CO-7/6 inch shsiffc 



frequency 


/' - *- — - — r> . « 1 1 

vJOnPtu 

‘*s>iHm& 


Coitgwied 


frequency 


Cbsereed 

\r.?liiud» 


fcccqxstedr 

ol itod® 


_ 










* 


200 


S.0 


4.2 


140 


2.5 


2.2 


230 


12.0 


13.3 


210 


9.5 


14.2 


240 


24.0 


23.5 


220 


54.0 


23.0 


245 


43.5 


41.0 


230 


62.0 


76-0 


240 


72.0 


76.0 


250 


20.6 


23.0 


235 


19.5 


31.0 


«. « 


11.0 


14.0 


275 


15.5 


18.7 


350 


5.0 


6.0 


300 


9.0 


10.0 








350 


5.5 


5.5 









•Jot©: "A" indicates operation on upper branch of non-linear response curve 
*j>" indicates operation cm lower branch of non-linear response curve 
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